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Abstract 



We study nonlinear ground states of the Gross-Pitaevskii equation in the space of one, 
two and three dimensions with a radially symmetric harmonic potential. The Thomas-Fermi 
approximation of ground states on various spatial scales was recently justified using variational 
methods. We justify here the Thomas-Fermi approximation on an uniform spatial scale using 
the Painlcve-II equation. In the space of one dimension, these results allow us to characterize 
the distribution of eigenvalues in the point spectrum of the Schrodinger operator associated 
with the nonlinear ground state. 

1 Introduction 

Recent experiments with Bose-Einstein condensates [PS] have stimulated new interest in the 
Gross-Pitaevskii equation with a harmonic potential. We take this equation in the form 

iu t + e 2 Au + (1 - \x\ 2 )u- \u\ 2 u = 0, x £ M. d , (eR + , (1.1) 

where the space dimension is d is one, two or three, u(x, t) £ C is the wave function of the 
repulsive Bose gas in the mean-field approximation, and e is a small parameter that corresponds 
to the Thomas-Fermi approximation of a nearly compact atomic cloud [Fer|, IT] . 

A ground state of the Bose-Einstein condensate is a positive, time-independent solution 
u(x, t) = rj £ (x) of the Gross-Pitaevskii equation (jl.ip . More precisely, r\ £ : W 1 \— > R satisfies 
the stationary Gross-Pitaevskii equation 

e 2 Ar] e (x) + (1 - \x\ 2 )t] £ (x) - rf £ (x) =0, x G R d , (1.2) 

i] £ (x) > for all x G R d , and r\ e has a finite energy E £ (r] e ), where E £ is given by 

E £ (u)= [ (e 2 \Vu\ 2 + {\x\ 2 -l)u 2 + \u A J dx. 

jR d V 2 / 

For d = 2, existence and uniqueness of a radially symmetric ground state r/ £ for a fixed, 
sufficiently small e > is proven in Theorem 2.1 of Ignat & Millot [IM| similarly to earlier 
works of Brezis & Oswald [BO] and Aftalion, Alama, & Bronsard [AAB] in bounded domains. 
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It is also shown in |IM] that rj e (x) converges to t]q(x) as e — > for all x £ M 2 , where r/o is the 
Thomas-Fermi's compact function 

, s / (l-|x| 2 )V2 for \x\ <1, n 

%(:E) = i far U| > 1. (L3) 

To be precise, Proposition 2.1 of [IM] states that for d = 2, e > sufficiently small, 

< i] £ (x) < Ce^exp f ^J^f ) for M > 1, (1-4) 

0^ (1 — |^| 2 ) 1/2 — 7? £ (^) ^ Ce 1/3 (1 - M 2 ) 1 / 2 for \x\^l-e 1/3 , (1.5) 



and 



1% - ?7o llci(iC) < C K e 2 , (1.6) 



where -fC is any compact subset of {x G M 2 : |x| < 1} and C and Ck are e- independent positive 
constants. The method used by Ignat & Millot in the case d = 2 to prove the existence of 
a radially symmetric ground state r/ E can be extended to the cases d = 1,3, even though the 
uniqueness of the ground state does not follow from |IMj for d = 3. We are concerned here 
with a uniform asymptotic approximation of the ground state rj £ on W 1 , in the limit e — > 0, for 
d= 1,2,3. 

At least two attempts have been made in physics literature |BTNN| KK to establish connec- 
tion between the nonlinear ground state rj £ f or d = 1 and solutions of the Painleve-II equation 

4i/" (y) + yi /(y) - ^ 3 (y) =0, (1.7) 

This equation arises as the formal limit as e — ► of the differential equation satisfied by v £ : 

4(1 - e 2 / 3 y yj(y) - 2e 2 ^dv' £ {y)+yv £ {y) - v 3 E (y) = 0, y £ (-oo,£- 2 / 3 ), 

where v £ is defined by 

rj £ (x) = e 1 / 3 u £ (y), y = 1 ~ 1 . (1.8) 

The convergence of ?y e to % as e -> suggests that we should consider the Hasting-McLeod 
solution uq of the Painleve-II equation |HM] , which is the unique solution of (jl.7p such that 

v o{y) ~ y 1 ^ 2 as V ~^ +°° an£ i ^o(y) ~~ * as y — > — oo. 

In both papers [BTNN, KK], the asymptotic solution ry £ is constructed at three spatial scales 

I : |x| ^ 1 — e 2 / 3 , II : \x\ G (1 -e 2/3 ,l + e 2/3 ), and III : |x| ^ 1 + e 2/3 . 

Solutions of the Painleve-II equation (|1.7|) are used at the intermediate scale II for matching 
conditions and connection formulas between the WKB solutions at the inner scale I and the Airy 
function solutions at the outer scale III. The same formal approach is also developed in [ZAKP] 
for approximations of excited states of the stationary Gross-Pitaevskii equation in the case d = 1 . 

We address the problem of uniform asymptotic approximations of the ground state r] £ of the 
stationary Gross-Pitaevskii equation (jl.2p using the Hasting-McLeod solution of the Painleve-II 
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equation (jl.7p . Our main result (Theorem [T]) in Section 2 establishes this approximation on a 
rigorous level. In the case when d = 1, we also study eigenvalues of the Schrodinger operator 

L% = -e 2 d 2 x + V £ (x), V £ (x) = 3r, 2 £ {x) -l + x 2 , 

that arises in the linearization of the stationary Gross-Pitaevskii equation (II. 2D at the ground 
state T] E . We prove in Section 3 that the spectrum of L £ + in L 2 (M) consists of an infinite sequence 
of positive eigenvalues {A^} n ^i such that for any fixed integer k ^ 1, 

Al fe ~ fi k e 2/3 as e^O, (1.9) 

where fi^ is the k th eigenvalue of the Schrodinger operator 

M = -4d 2 + W (y), W (y) = 3u (y)-y. 

We note that Mq arises in the linearization of the Painleve-II equation (|1.T[) at the Hasting- 
McLeod solution vq. Therefore, the scaling transformation (|1.8j) leading to the Painleve-II equa- 
tion (jl.7p becomes useful for analysis of eigenvalues of the Schrodinger operator IF,. 

It is clear from the shape of r\ e that the operator L £ + has a double- well potential V e (x) with two 
symmetric minima converging to ±1 as e — > 0, while the operator Mq has a single- well potential 
Wo(y). These facts explain both the asymptotic correspondence between eigenvalues of L e + and 
Mq and the double degeneracy of each pair of eigenvalues in the asymptotic limit (|1.9|) . Formal 
results of the semi-classical theory for the operator L e + are collected in Section 4. 

While a different technique is exploited in our previous work [GPj . the result (jl.9p provides 
the same kind of asymptotic behaviour for the smallest eigenvalue of L £ + as the one we obtained 
for the lowest eigenvalue of the simplified operator 

L% = -e 2 d 2 + V (x), V {x) = 3 V 2 {x) -l + x 2 . 

The spectral stability of the ground state in the Gross-Pitaevskii equation (jl.ip is deducted 
from the analysis of the symplectically coupled eigenvalue problem for Schrodinger operators L e + 
and LI, where 

LI = -e 2 dl + V £ (x), V e {x) = ri 2 (x) - 1 + x 2 = 

Unfortunately, the asymptotic scaling (|1.8p leading to the Painleve-II equation (|1.7p does not give 
a correct scaling of the eigenvalues of L e _ nor the eigenvalues of the spectral stability problem 
because the potential V £ (x) is a single well with a nearly flat bottom on the interval [—1,1], 
which is mapped to [0, e~ 2 ^} by the change of variable y = (1 — x 2 )/e 2 ^. Analysis of the 
eigenvalues of the spectral stability problem and construction of excited states of the stationary 
Gross-Pitaevskii equation are two open problems beyond the scope of this article. 

Notations. If A and B are two quantities depending on a parameter e belonging to a neigh- 
borhood £ of 0, 

• A(e) < B(e) indicates that there exists a positive constant C such that 

A(e) ^ CB{e) for every e e £. 
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• A(e) ~ 5(e) if A(e)/B(e) -> 1 as e -> 

• A(e) = 0{B(e)) as e —> if A(e)/B(e) remains bounded as e — > 0. 

Let F{x) be a function denned in a neighborhood of oo. Given a G R, {/ m } m gN G R, and 
7 > 0, the notation 

.r — >no * — • 

m=0 

means that for every M G N, 

M 

F(x) - x Q ^ / m x" 7m = 0(x Q -^( M+1 )) as x -» oo, 

m=0 

and, moreover, that the asymptotic series can be differentiated term by term. 
We use the following spaces: 

• H°°(R) = n H S (R), where H S (R) is the standard Sobolev space. 

• L^(R d ) is the subspace of radially symmetric functions in L 2 (R d ). Note that if /(| • |) G 
L 2 (M. d ), then 

/>oo 

3<2-li / „d-l| r c„m2. 



ll/(H)lb (K «) = F^l / r*- 1 !/^)!^, 





where |S rf_1 | is the surface of the unit sphere in R d . Similarly, \M d \ is the volume of the unit 
ball in R d . 



2 Uniform asymptotic expansion of r) e 

In what follows, d = 1, 2 or 3 and e > is sufficiently small such that, as it is proved in Theorem 
2.1 of |IM] , there exists a positive classical solution r] £ of 

e 2 A?] e (x) + (1 - \x\ 2 )t] £ (x) - r)l(x) =0, x G M d . (2.1) 

Moreover, this ground state % is radially symmetric, so that we can define a function ^ e on 

J £ := (-oo,e- 2 / 3 ] by 

r ?£ (x)=e 1 / 3 ^(^j^), xGR d . (2.2) 

Let y = (1 - \x\ 2 )/e 2/?J be a new variable. Notice that y covers once J £ as \x\ covers R+. It is 
equivalent for r] E to solve (|2.ip and for z^ e to solve the differential equation 

4(1 - £ 2/3 2/K'(y) - 2e 2 / z du' £ {y) + yu e (y) - v 3 £ (y) =0, y G J e . (2.3) 

Let TV ^ be an integer. We look for v £ using the form 

N 

u £ (y) = ^2e 2n ^ n (y)+e 2 ^/ 3 R N ,M, y G J £ . (2.4) 

n=0 

Expansion (|2.4|) provides a solution of equation (|2.3|) if {z4i}o^n<7V and Rn,e satisfy equations 
(1231) and <|2ZD below. 
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uq solves the Painleve-II equation 

H(y) + vMy)-vl(y) = ^ ^r, 

for 1 ^ n ^ N, v n solves 

-K(y) + Wo(y)v n (y) = Fn{y), i/eR, 

where 



(2.5) 



(2.6) 



and 



Rn e solves 



W (y) = 3^(y) - y 



F n(y) = ~ "nAy^miy^miy) - 2dv' n _ 1 [y) - Ayv'^y), 



n-i + n 2 + n 3 — n 



4(1 - e 2 ' z y)R" N£ + 2s 2 'HR' n £ + W R N>S = F N)£ (y, R N)£ ), y G J e , (2.7) 



where 



27V- 1 



F Nte (y,R) = -(4y^ + 2du' N )- £ e 2 ™/ 3 ]T 



2/V 



\ 



71=1 n-i + ri2 — n 



"2 



5^ 72. 1 , 712, 71 3 "S: 

2N+1 



R-ls e^ n _ (N+l) \R 2 -e^^R\ 



J 



=N+1 



Notice that for ^ n ^ iV, ^ n (y) is defined for all y E R and does not depend on e, whereas 
RN,e(y) is a priori only defined for y £ J £ . 

Appropriate solutions of system (|2.5|) , (|2.6p and (|2.7|) enable us to prove the following theorem. 



Theorem 1 Let vq be the unique solution of the Painleve II equation \2. 5\) such that 

v o{y) ~ y 1 ^ 2 as y ~~ * +°° an d u o{y) ~^ o as 2/ ~ *• —00. 

For n ^ 1, i/iere exists a unique solution v n of equation \2. 6\) in H°°(M). For every N ^ 0, there 
exists en > and Cat > such that for every < e < en , there is a solution Rn,e £ C 00 nL 00 (J e ) 
o/ equation {2. 7\ ) with 

< Civ^ ( ^ 1)/3 and x ^ P'jgH G // 2 ' ' 



A' 



2n/3 



n=0 



1 — |x| 

^27F 



1 — \x\ 



, x G 



is a ground state of equation \2.1\ ). 
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Remark 2.1 For d = 3, the remainder term in \2. 8\) may have the same order as the last term 
in the sum, because of the growth of the upper bound on ||-Riv,e 1 1 z,°°(j e ) as e I 0. 

Remark 2.2 For d= 1,2, the ground state we find in Theorem{l\is the unique ground state of 
equation ( fff.ij) . thanks to the uniqueness result proved in lIMjj . For d = 3, it is not clear whether 
the ground state of Theorem^ coincides with the one obtained by the method of Ignat and Millot 
in JIMj . because uniqueness of a ground state does not follow from JIM) /. 

The proof of Theorem Q] is described in the following three subsections. Notice first that it is 
sufficient to prove the Theorem for an arbitrarily large value of N. Indeed, for every integer 
Nq > 0, the result of the Theorem for N < Nq is a direct consequence of the result for N = Nq. 
Also, for convenience, we shall assume in the sequel that N ^ 2. 



2.1 Construction of v n for ^ n ^ iV 

We are looking for a solution v e {y) of equation f|2.3|) that satisfies the following limit as e —* 0: 

E V V(£ -/3 (1 _ l2)) _( U-;f) 1/2 ^ xg[-1,1], 
v v " e^o \ for \x\ ^ 1. 

Therefore, we choose vo(y) to be the unique solution of the Painleve-II equation (|2.5p that satisfies 
the asymptotic behavior v$(y) ~ y 1 / 2 as y — > +oo and converges to zero as y — * — oo. Existence 
and uniqueness of this solution are proved by Hastings & McLeod |HM| . Asymptotic behaviour 
of vo(y) as y ^ ±oo is described in more details in Theorem 11.7 of [FIKN]. These results are 
combined together in the following proposition. 

Proposition 2.1 ]HM\ \FIKfflj The Painleve-II equation 

Au"(y)+yu{y)-u 3 (y)=0, y E R, 

admits a unique solution vq G C°°(R) such that 

u o(y) ~ y 1 ^ 2 a s y —> +°o and vo(y) — ► as y — > — oo. 

Moreover, vq is strictly increasing on R ; v'q has exactly one zero on R, which is an inflection 
point of vq. The behaviour of vq as y — > — oo is described by 

, ( y ) = -^(-2 y )- 1 /Vl(-^ 3/2 (l + O(| y |- 3 / 4 )) » 0, (2.9) 
2y7T V / y^-oo 

whereas as y — > +oo, it is described by 



oo , 

On 



where bo = 1, b\ = 0, and /or n ^ ; 

^ n+1 n+1 n+2-Z 

K+2 = 4(9n 2 - l)6 n - - ^ bmbn+2-m ~ ^ X) X] h lb m b n + 2 -l-n 

m=l 1=1 m=l 
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Next, we construct v n G H°°(R) for n ^ 1 by induction on n. For n ^ 0, we consider the 
following property: 

f • v k G F°°(R) solves ([22]) (with n replaced by fc), 

VkW) ~ /~ 2fc Em=oa,™2/~ 3m/2 for some {ft.mlmeN, 

^fc(y) ~ o, 

2/->-oo 



(H n ) Vfce {!,..., n}, ^ 



where 

-5/2 if d = 1, 
1/2 if d = 2, 3. 

(Ho) is empty and, therefore, true by convention. Fix n 1 and assume that (H n _i) is true. 
We are going to construct v n such that (H n ) is satisfied. We will make use of the following two 
lemmas, which are proved in Sections 5 and 6. 

Lemma 2.1 Let W G C X (R) such that W € L°°(IR + ) and there exists C ,C + , A + > with 

W(x) ^ C + x for x ^ A+, W(x) > C for x e R, and W'{x) ^ for x ^ A+. 
lei / G L 2 (R) such that x a f G L°°(A + , +oo) /or some a > 0. Let 

ip = (-d 2 x + Wy 1 f£H 1 (R). 

Then, as x — > +oo, 

ip{x) = 0{x- {a+1) ). (2.11) 
Moreover, if f and W admit asymptotic series 

+oo +oo 

f(x) « x- Q Vc m x^ m , W(x) « xV Um x^ m , (2.12) 

x— ++oo • 2:— >+oo ' 

m=0 m=0 

/or some coefficients {c m } mg N, {%}meN Qftd 7 > such that 3/7 is an integer, then (p admits 
an asymptotic series 

+00 

tp(x) « V d m x~ 7 "\ (2.13) 

:r— >+oo ^— — ' 

m=0 

/or some coefficients {d m } me N- iVi particular, as x — > +00 ; 

y?'(x) = 0(x-( Q+2 )), </(x) = 0(x"( a+3 )). (2.14) 

Lemma 2.2 Zei Wo(y) := 3fg(y) — y, where vo(y) is the solution of the Painleve-II equation 
\2. 5\ ) given in Proposition ^. 1\ Then, 

W min := inf W (y) > 0. 
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^,From the asymptotic behaviors of vo(,y ) as y — > ±00, we infer that 

Wo(y) ~ 2y as y — > +00 and Wo(y) ~ — y as y — ► —00. (2-15) 
Let us consider the operator 

M :=-Ad 2 y + W (y) 

on L 2 (R) with the domain, 

Dom(Mo) = {u G L 2 (R) : -4«" + W u € L 2 (R)}. 

The Schrodinger operator Mq arises in the linearization of the Painleve-II equation at v = uq. 
The spectrum of Mq is purely discrete and, thanks to Lemma 12.21 it consists of a sequence of 
strictly positive eigenvalues which goes to infinity. If n = 1, it follows from the choice of vq and 
from properties (jlH)]) - (f2TT0|) that 

+00 

F x {y) « (1 - d)<T 1/2 + y- 7/2 V 3(m + 2)6 m+2 2- 3 ( m+2 )/ 2 (l - 3(m + 2))y~ 3m / 2 (2.16) 

7/ — >+oo ^ * 



y— >+oo 

m=0 



iq(y) « 0. (2.17) 



and 



Thanks to Lemma 12.2^ we can look for z^i solution of (|2.6p with n = 1 in the form 

(1 - d)<Z>(y) 

where <£ G C°°(R) is such that 3>(y) = if y ^ 1/2, $(y) = 1 if y ^ 1. Then, v\ has to solve 

- 4u'{{,y)+W (y)h(y) = h{y) ■■= F x { y ) - (1 - d)y- l ' 2 <s>(y)+^ ( ( ^~ff/2 } ) > f G M -( 2 - 18 ) 

^From the asymptotic expansions (|2.9p - (|2.10p of i/o, we infer that Wo also admits asymptotic 
expansions as y — ► ±00. Since moreover vq G C°°(R), it follows from (12.151) . (12.16p . (12.17P and 
(|2.18p that F\ G i5f°°(R). Then, property (Hi) follows from Lemma |2 . 1 1 applied on the one side 
to v\ := M Q ~ 1 Fi with a = 7/2, so that v\(y) = 0(y~ 9 / 2 ) as y — > +00, and on the other side to 
y l— ► vi(—y) with a arbitrarily large. Furthermore, if n ^ 2, we have 

-^n(y) = - l/ ni(y)^n 2 (y)^n 3 (y) - 3 ^ v (y)v ni (y)v n2 (y) -2dv' n _ 1 (y) - Ayv'^iy). 



< 711,712,713 < 77 < 771,772 < 77 

771 + 77 2 + 773— 77 771 + 772 — 77 



Thanks to (H n _i), all the terms in the right hand side admit an asymptotic expansion at ±00. 
More precisely, 

+00 

F n (y) « y p+1 ~ 2n £ fn, m y~ 3m/2 

y— >+oo » 
m=0 



for some coefficients {f n ,m}meN, whereas 



F n (y) ~ 0. 
2/-+-00 
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Since F n G C°°(M.) and n ^ 2, we deduce that F n G J?°°(R), and we can define v n = M Q 1 F n G 
fP°(R). By Lemma O with 7 = 3/2 and a = 2n - /3 - 1, we then have 

+00 

m=0 

for some coefficients {g n ,m}meN, and 

where we have applied Lemma, 12. 1 1 with, / y — 3/2 to the function z^ n ( — ?/). Therefore, (H n ) is true, 
which completes the construction by induction of the sequence of solutions {v n {y)} n ^\ of the 
inhomogeneous equations (j2.6j) . 



2.2 Construction of Rm )£ 

In this subsection, we construct a solution Rn,e to equation (|2.7p . such that given the ^ n 's 
constructed in subsection 2.1, expansion (|2.4jl provides a solution of equation (12. 3p . The solution 
i?jv,e of equation (|2.7p is obtained by a fixed point argument. In order to explain the functional 
framework in which the fixed point theorem will be applied, let us first introduce the functional 

spaces 

L 2 £ = {u G L\ oc (J e ) : (1 - e^yfl^u G L 2 (J £ )} 

and 

H] = {u G L 2 : (1 - e 2 ^y) dl ^ G L 2 (J £ ) and (1 - e^yf/^W^u G L 2 (J e )} , 
endowed with their respective squared norms 



re" 2 / 3 

|2 



lull- := Hi, = (1 - e 2 / V /2 " Vdy 



and 2 a 

[4(1 - e 2 ^y) d ' 2 \u'\ 2 + (1 - ^2/3^^/2-1^^21 

') 

We are looking for a solution R^ y£ {y) of Equation (|2.7p on J £ such that the function R^^{e~ 2 ^{l- 
\x\ 2 )) is regular on M. d . As a result, it is convenient for the sequel to introduce the map T £ : L £ \- 
L 2 {R d ) defined for u G L 2 £ by 

(T £ u)(z) :=u(e- 2 / 3 -e 2 / 3 \z\ 2 ), 



which makes the link between functions defined on J £ and radial functions defined on M, d , in 
terms of the variable z = e~ 2 ^x G M d . An easy calculation shows that T £ is a bijection from L 2 
into L 2 (M. d ), and that for every u G L 2 , 



2 l S<i_1 | 2 

|T £ u|| i2(Rd) = 2e 2(d-l)/ 3 IMIe- ( 2 - 19 ) 
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Moreover, T £ induces a bijection from into 



Q £ := \u G L 2 (R d ) 
and for every u G H £ , 



\Vu\ 2 + W (e- 



-2/3_ £ 2/3| z |2 )M 2 



\VT £ u\ 2 + W (e- 2/3 - e 2 / 3 \z\ 2 )\T £ u\ 2 \ dz = ^^J uf^. (2.20) 



dz < oo 



Let us rewrite equation (|2.7p for the remainder term Rn >£ (jj) in the operator form 

M £ R N , £ (y) = F N , £ {y, R NjE ), y G J £ , (2.21) 
where M e is the self-adjoint operator on L 2 defined by 



M £ := -4(1 - e 2 / 3 y)- d / 2+1 d y (l - e 2 / 3 y) d / 2 d y + W (y) = {T e )~ l K e T e 
Dom(M e ) = {u G L 2 : K £ T e u G L 2 (R d )} 

and denotes the Schrodinger operator on L 2 (M. d ), 

K e := _ A + ^ o(e -2/3_ e 2/3| z | 2)> 



(2.22) 



The solution -R/v,e of the nonlinear equation (|2.2ip will be obtained from the fixed point theorem 
applied to the map 

&N,e '■ R ^ {M £ )~ 1 FN, e {-, R), 

which will be shown to be continuous from H £ into itself. First, we shall prove the following 
lemma. 

Lemma 2.3 The operator M £ is invertible, and for every f G L 2 , 

IKMTVlk ^"ifll/ll,. 

Proof. Let us consider the continuous, bilinear, coercive form on Q £ defined by 



a(u, v) 



VuVv + W {e 



-2/3_ £ 2/3| z |2 )TO 



dz. 



By the Cauchy-Schwarz inequality and Lemma 12.21 for every / G L 2 



J 



T £ fvdz 



defines a continuous linear form on Q £ . Thus, by the Lax-Milgram Theorem [GTj . there exists a 
unique tp G Q £ such that for every v G Q £ , 



a(ip,v) 

Moreover, ifi G Q £ is radial and satisfies 



T £ fvdz. 



K £ tp = T £ f in V 



i tiad\ 
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Thus, <p := (T 6 )- 1 ^ e Hi n Dom(M £ ) satisfies 

M £ ip = f. 

^From (|2.20p and a calculation similar to (|2.19p . we also check that 



2e 2(d-i)/3 2e 2 ( d - 1 )/ 3 



-2/3 



-1/2 



elMIe < Wmin H/lkllVllffl- 



from which the upper bound on tp = (M £ ) /in Hi follows. 

Next, we prove that R i— ► Fjy j£ (-,i2) continuously maps if* into L 2 . We write 

ify e (y, fl) = F Nfi {y) + G N , e (y, R), 

where 



F N ,o = -{±yv% + 2dv' N )- Y, 



(2.23) 



n\ + n 2 + n 3 — AT + 1 
^ Tii j 71 2 ! 71 3 ^ 



and 



G 



2N-1 



2iV 



V 



3j> 2n/3 E 

n=l m + n2 — n 

< ni , n 2 ^ N 



- E e2n/3 E 

71 = 1 ni + n 2 + ri 3 = n + A 7 + 1 

^ ni , no , r? 3 $J N 

( 2N+1 \ 

- 3 E ^ /3 -n-(^ + l)U 2 -^ +1)/3 ^. 
\ n=jV+l / 

We first show that .F/v,o £ -^e- Indeed, from the properties of the u n 's, we infer that 

F N ,o(y) ~ o. 

and i^jv.o also admits an asymptotic expansion as y — > +oo, with 

Ayp'M + 2dp' N {y) = 0{yP- 1 ~ 2N ) 

and if ni + tt-2 + r^3 = A r + 1, 



i? 



(2.24) 



v ril v ri2 v nz (y) = < 



' 0(y- 19 / 2 - 27V ) if d = 1 and m, n 2 ,n 3 > 0, 
C(y ~ 13 / 2 ~ 2 -' v ) if d = 1 and ni or 712 or 71,3 = 0, 

{ 0{y -l/2-2N } if02) 



(notice that n\ + ri2 + = N + 1 with ^ Hi,rt2, rt3 ^ A r implies that at most one of the 
numbers ni,ri2,Ti3 is equal to 0). Since N ^ 2, we deduce that in any case, 

F Nfi (y) = 0{y- g / 2 ) as y - +00, while ityofa) ~ 0. 
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Therefore, for a > sufficiently large and e < 1, 

-2/3 i 



f (l-e^yY^F^dy < [ (l + \y\T 2a dy+ f ~ y~\l - e^yf^dy. 

J —oo J —oo Jl 

In the case d = 1, the second integral in the right hand side is estimated by 



y "(1 ~ 2:1 //* l/2 dy < £ 16/3 / 77^17^ 

7 e 2/3 (l - zyi 1 



l z -9 



rl/2 -16/3 rl i 



whereas for d ^ 2, 

— 2/3 — 2/3 

2/- 9 (l - fvf^dy < jf" y- 9 dy<l. 
Therefore in both cases Fjv,o £ L 2 and 

ll^olle < 1- (2.25) 

Similarly, the term which does not depend on R in the right hand side of (|2.24p is 0^2 (e 2 / 3 ). 
Let R G . To estimate the linear term in R in the definition of Gjv,e, notice that if n\ + ri2 = 
n ^ 1, then ni or 112 is not equal to 0, thus v ni v n2 {y) = 0{y~ l ) as y — > +oo. In particular, 
v ni Vn 2 S and 

— 1/2 

In order to estimate the quadratic and cubic terms in the right hand side of (|2.24p . the following 
lemma will be useful. 

Lemma 2.4 Let p = 1, 2 or 3. There exists a e -independent constant C > such that for every 
e > 0, if u £ Hi then u v £ L? e and 

£ 

Proof. Let ueH^. We have checked in (ET201) that T e u eQ £ C H l (R d ) and 

l|T e u||^ 1(Rd) < ||r £ n|| Qe < e-W 3 \\u\\ H i. (2.27) 

By Sobolev embeddings, it follows that T £ u £ L 2p (M. d ), and 

|| U P|| £ < £ ( d - 1 )/ 3 ||T e K)|| i2(Kd) =e (d - 1)/3 ||T £ n||^ (Rd) 

< £ (^D/3|| T£u ||^ (Rd) < e -(P-l)(^l)/3|| u ||^, ( 2 . 28 ) 

where we have also made use of (|2.19p with u replaced by vP . m 
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Remark 2.3 The statement of Lemma \2.J\ can be extended for all values of p for which i/ 1 ! 
is continuously embedded into L 2p (M d ), that is 1 ^ p ^ oo for d = 1, 1 ^ p < oo for d = 2 and 
l^p^3ford = 3. 

Thanks to Lemma 12.41 f° r an Y integer k ^ 1, 

HR% < ll^ll^^e-^ 1 )/ 3 !^!!^ < e-^-^WRW 2 ^, (2.29) 

whereas for k = 0, 

HR% < |ko||Loo(j £ )£- { ^ 1)/3 ||ii||ii < £- d/3 ||i?|lii (2.30) 

On the other side, 

\\R 3 \\s<£~ 2{d ~ 1)/3 \\R\\Hi, (2-31) 

thanks to Lemma E3 again. By Lemma E31 as well as bounds (f2T25j) . (|2T26j) . ([239|) . ([230]) and 
(I2T3TH . 

N^^v^C^) - < £ 2/3 + ^ 2 / 3 ||^||^ + ^C 2iV + 2 - rf >/ 3 || J R||^x ^ ^C 4iV + 6 - 2rf >/ 3 ||J^||^ 1 , 

where 

R% >£ : = (Mr 1 ^. 

In particular, for e > sufficiently small and for some e-independent constant C > 0, &n,s maps 
the ball 

B e :=B Hl {R%^Ce 2 ^) 

into itself, where we have used the assumption N ^ 2. Similarly, there exists an e-independent 
constant C > such that for every R\, i?2 in B £ , 

H^ORi) - $ N> e( R 2)\\Hi ^ Ce^WRr - R 2 \\ H i. 

As a result, provided e is sufficiently small, &n,e is a contraction on B e . The Fixed Point Theorem 
ensures that &N,e has a unique fixed point Rjy,e £ B £ . In particular, 

\\RN,e-R N jHl <e 2/3 . (2.32) 

We next prove that Rn, £ satisfies the regularity properties stated in Theorem [TJ The fixed point 
Rn,e £ He °f ^N,e has been constructed in such a way that T e R^ )£ G H 1 (M d ) solves the equation 

K E T e R Ni£ = T%F N>£ (-,R N>£ )) E L 2 (R d ). (2.33) 

Thanks to Lemma 12.31 and (I2.25|) . we obtain 

lK, £ lki<ll^o|| £ <l. (2.34) 

Thus, (I2T32|) yields 

\\RN,e\\e<\\RN,e\\m<l- (2-35) 
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As a result, from (|2^5|) . (12T26J) . (l239j) . (|230|) and ([23TD . we infer 

\\F Nj£ (-,R Nj£ )\\ £ <1. (2.36) 
^From (12331) . (12361) and (12T91) we deduce 

ll/rr^H^) = iit^f^^.,^))!!^^) < e-^ 3 . (2.37) 

Next, we use the following Lemma, which is proved in Section 7. 
Lemma 2.5 (K 6 )^ 1 G £(L 2 (M d ), # 2 (M d )) is uniformly bounded in e. 

As a result, we infer from the Sobolev embedding of H 2 (R d ) into L°°(lR d ) that l^i^e G H 2 (R d ) 
and 

||-RAr,£|U°°(j e ) = ||T e i?iv,e||ioo( R d) < e _(d ~ 1)/3 . (2.38) 
Moreover, a bootstrapping argument shows that T £ Rn >£ G C°°(1R c( ). As a result, Rn,e £ C°°(J e ). 

2.3 i/ s (j/) > for all y e J e 

We have constructed above {v n }n^o and Rn,s i n such a way that 

N /I UI2' 



/ 1 I 1 2 \ 

£ l/3^ £2 n/3^ (±ZM-j +^/ 3 + 1 (T^^ £ )( £ - 2 / 3 x) 



T] e {X 

n=0 

= e 1/3 ^ e (^ajr-) . ( 2 -39) 

is a classical, radially symmetric solution of equation (|2.ip . In order to claim that fj £ is a ground 
state, it is sufficient to check that fj £ (x) > 0, for every x G M d , which is equivalent to v e {y) > 0, 
for every y G J £ . 

For every n ^ 1, Hi'nlli^fR) i$ 1- Therefore, from (|2.38j) . (|2.39j) . since N ^ 2, we deduce the 
existence of a constant C > such that for every y G J e , 

f e (y) - ^o(y) ^ -Ce 2 / 3 . 

Since fo(y) increases from to +oo as y goes from — oo to +oo, we deduce that for e <C 1, 

u e (y) > ^o(-l) - Ce 2 ' 3 > 0, y G [-l,e" 2/3 ]. 

Coming back to the variable x, it follows that 

fj £ (x) > 0, |a?| < (1 + e 2/3 ) 1/2 . (2.40) 

It remains to prove that %(x) > for all |x| > (1 + e 2 / 3 ) 1 / 2 . Assume by contradiction that fj £ is 
not strictly positive on M. d . Then, let 



inf{r > 0, fj £ (r) = 0} G ((1 + e 2/3 ) 1/2 , oo) 
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where for convenience, since fj £ is radial, we denote = fj £ {x). By construction, fj £ {r £ ) = 

and fj' £ (r £ ) ^ 0. If fj' £ (r £ ) = 0, then fj £ = 0, because r? e (r) satisfies the differential equation 

-^iTr y& e ) (r) + h {r " " 1 + ^(r) 2 )^) = °- 
This is a contradiction with (|2,40p . Thus, fj' £ (r £ ) < 0. Let 

f £ := sup{r > r £ ,fj £ (r') < for r' G (r e ,r)} G (r e ,+oo]. 
Then, for every r G (r £ ,f e ), 

A (V" 1 ^ ( r ) = ^(r 2 - 1 + Ve(r) 2 )Ur) < 0, 
ar \ ar J e A 

and we deduce by integration that for every r G (r e ,f E ), 

r d -^(r) <rf"V e (r e ) < 0, 

and 

fj £ (r) ^ r d £ - l i)' £ (r £ ) f s^ds. (2.41) 

The right hand side in (|2.4ip is a negative, decreasing function of r, which implies r £ = +oo, as 
well as a contradiction with the fact that fj £ (r) — » as r — > +oo. Therefore %(r) > for all 
r G M+. 



3 Spectrum of the Schrodinger operator L £ + in the case d = 1 

Consider the Schrodinger operator 

= ~e 2 d 2 x + V £ (x), V £ (x) = to£(x) - 1 + x 2 , 

associated with the stationary Gross-Pitaevskii equation (jl.2p linearized at the ground state rj £ . 
It is a self-adjoint operator on L 2 (R). Since the potential V £ (x) is confining in the sense of 
V £ {x) — > +oo as \x\ — * oo, L £ + has compact resolvent and a purely discrete spectrum. By Sturm- 
Liouville theory, the eigenvalues of IA, denoted {A^} n j>i (sorted in increasing order) are simple. 
Moreover, thanks to the even symmetry of V £ on R, the eigenfunctions of L £ + corresponding to 
A^j are even (resp. odd) in x if n is odd (resp. even). If A is an eigenvalue of L £ + and (p G L 2 (R) 
is a corresponding eigenfunction, we define a function v G L 2 by 

<P(z) = V ^2/3 1 , x G M+. 

Let us denote W e (y) = 3v 2 (y) — y. Then, <£> G L 2 (M) is an even eigenfunction of L £ + corresponding 
to the eigenvalue A if and only if v G L 2 satisfies the differential equation 

(_ 4 (i - e^yf/ 2 d v {l - e^y^dy + W £ {y)) v(y) = e^My), y G J £ (3.1) 
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and the Neumann boundary condition 

c//(0) = -2e~ 2 / 3 ((1 - e^yY^v'iy)) = 0. (NC) 

y=e~ 2 / 3 

Similarly, tp £ L 2 (M) is an odd eigenfunction of L £ + corresponding to the eigenvalue A if and only 
if v £ L 2 satisfies (|3.ip and the Dirichlet boundary condition 

y,(0) = v(e- 2 / 3 ) = 0. (DC) 
As a result, the eigenvalues of L £ + are directly related to the eigenvalues of the two self-adjoint 
operators on L 2 , M £ and M £ , where 

M £ := -4(1 - e 2 / ?J y) l / 2 d y {\ - e 2 l 3 y) l l 2 d y + W £ {y), 
Dom(M e ) = {v £ L 2 e : M £ v G L 2 and v satisfies (NC)} , 

and M £ is defined similarly by replacing (NC) by (DC) in the definition of the domain. Namely, 
if we denote {Anln^i (resp. {/^}n^i) the eigenvalues of M £ (resp. M £ ) sorted in increasing 
order, then for every n ^ 1, 

An = £~ 2/3 >%n-l and /in = ^ 2/3 A!„. 

As s — > 0, the eigenvalue problems (|3.ip for the operators M £ and M £ formally converge to the 
eigenvalue problem for the Schrodinger operator Mq defined after Lemma 12.21 

(-4d 2 + W Q {y))v(y) = fiv(y), yeR, where fi = e~ 2/3 A. 

By the discussion below Lemma 12.21 the purely discrete spectrum of Mq in L 2 (R) consists of an 
increasing sequence of positive eigenvalues {/i n }n^i- We shall prove that the eigenvalues of L £ + 
converge to the eigenvalues of Mq as e — > 0, according to the following result. 

Theorem 2 The spectrum of L £ + consists of an increasing sequence of positive eigenvalues {A^} n ^i 
such that for each n ^ 1, 

lim%i=lim% =/ u n . (3.2) 

Proof. We prove only the convergence of fi £ n = A| n /e 2//3 to fi n , for every n ^ 1. The proof of 
the convergence of fi £ n = A|„„i/e 2 ^ 3 to /x n is identical. 

Denote by (•, •) and || • || the scalar product and the norm in L 2 (]R), and by (•, •) and || • || £ 
the scalar product and the norm in L 2 . lfu,ve L 2 (R), u _L v means that (u, v) = 0, whereas if 
u,v £ L 2 , u _L e v means that (u, v) £ = 0. We denote by 



the Rayleigh quotient for the operator M £ , where Q £ denotes the corresponding bilinear form 

Q £ (u,v) = jf (4(1 - e^y^dyudyv + _^ )1/2 «foMy)) e&, 
defined for u,v *E H^. Similarly, 



fl W = II..M2 
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denotes the Rayleigh quotient for Mq, where Q is the corresponding bilinear form 



Q(u,v) 



(AdyudyV + W (y)u(y)v(y)) dy, 



defined for u,v G {u G H l (R) : Wq /2 u G L 2 (M)}. 

Let u n (resp. u n ) denote an eigenfunction of M £ (resp. Mq) corresponding to the eigenvalue 
fi £ n (resp. fi n ), normalized by ||tt^|| e = 1 (resp. \\u n \\ = l).The eigenvalues of Mq are given by the 
Max-Min principle: 



inf R(v), 
v G Dom(Mo) 

V J_ Ui,- ■ ■ ,Un-i 



(3.3) 



whereas the eigenvalues of M £ are similarly given by 



inf 

v G Dom(M £ ) 
v ± E «!,•■■ 



(3.4) 



Let us fix (5 G (0,2/3). Let $ G C°°(M) be an non-decreasing function such that $ = on M_ 
and $ e 1 on [1, +oo). For e > sufficiently small, we also define Xe G C£°(M) by 



2x + e 



-2/3 



e-2/3 _ 2e - 



/ g- 2 / 3 - 2x 
£ -2/3 _ 2e -i 



-2/3 p-2/3. 



such that Xe is even, = 1 on [— e , e ] and Supp(x £ ) C [- , . 2 
recursively the following properties: 

(i) n ji n = fi n + 0{e^- 5 ), 

(ii) n for every fc > n + 1, (x £ Uk,u n ) £ = ©(e 1 / 3 ^/ 2 ), 

(iii) n for every fe > ra, (x e «|, u„._i) = ©(e 1 / 3 " 5 / 2 ), 

(iv) „ inf ||x £ n £ - cu„|| = ©(e 1 / 3 " 5 / 2 ), 

(v) n inf || Xe «n-l " CU n _ x \\ £ = 0(£ 1 /3-5/2) 5 



We shall prove 



(Gn) 



where for n = 1, (iii)i and (v)i have to be understood as empty properties. Let us fix n ^ 1 
and assume that (Gk) is true for every k G {1, • • • , n — 1} (for n = 1, this condition is empty, 
therefore true by convention). The proof of (G n ) is then divided in five steps. 



Step 1. Upper bound on fi e n . First, we shall prove that 



R E (V £ r , 



fi, 



+ 0{e 2/3 ~ s ), where 



n-l 



XeU r , 



^2 (^n, 



(3.5) 



k=l 



Then, thanks to (|3.4p . since v n G Span(uf , • • • } u £ n _ l )- Le C L 2 by construction, (|3.5p yields 



fi n ^ fi n + 0(e 



2/3-5^ 



(3.6) 
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^From (i)^ and (ii)k, which are satisfied for k ^ n — 1 thanks to the recursion assumption, we 
have 

n-l 2 

ytx..^)- Eft yfcl ,„ ae „„ ) + 0(e ,/3- a) 

* w = — — = llx^lll + ot^) ' ( ' 

k=l 



Next, 



-2/3 

,2„,2 



I ||2 _ / 2 Xg^n , 

IX^n|| £ - J g _ 2/3 (1 _ e y 3y y /2 a >y 



(1 + 0(62/8-')) T n 2 dy+ | i/ 3 "u/2 ^- ( 3 - 8 ) 



/ 2 /a TTWW* < ^ / u 2 dy < exp(-2e- 5 ) < e 2 / 3 , (3.9) 



The last term in the right hand side of (|3.8p is estimated as follows 

le "n 
?Zi (1 - e 2 / 3 y) ] 

where we have used the following Lemma. 

Lemma 3.1 For every m ^ 1, there exists a constant C m > such that for every y G R, 

|t*m(y)| < C m exp(-|y|) (3.10) 

and 

K(y)KC m (|y| + l)exp(-|y|). (3.11) 

Proof. Since Wo(y) — > +oo as y — ► oo, we can fix b n > such that inf{Wo(y) : |y| ^ & n } > 
4 + fi n . Then, 

(-4S 2 + W (y) - fi n )e-\y\ = (W {y) -fi n - 4) e -M > 0, |y| > 6 n . 
Since u n solves the eigenvalue problem 

(-4d 2 + W (y) - fi n )u n = 0, y G R, 
thanks to Corollary 2.8 in PQ, there exists C > such that 

K(y)K<? e -w, \y\^b n + l. 

Bound (IBTTUI) follows, since u n G Dom(M ) C i? x (R) C L°°(R). Then, from the differential 
equation M$u n = fi n u n and thanks to the asymptotic behaviour of Wq, we infer 

\<{y)\ = \Wu n {y)-W Q (y)u n {y)\<{\y\ + l)e-\v\, y G R. (3.12) 
By integration of (|3.12p between — oo and y, we deduce, for y < 0, 



KG/) I 



u n (s)ds 

-oo 



<(M + iK 



The same kind of estimate is obtained for y > by integration of (13.120 between y and +oo, 
which provides (|3.1ip and completes the proof of Lemma 13.11 ■ 



18 



Using Lemma 13. II again, as well as the normalization of u n , we infer that 



u 2 n dy = l + 0(e 2 ' 3 ). 

-s 

iFrom (pT8j) . (|33j) and (JHH]), we deduce that 

\\x £ u n f £ = l + 0{e 2 /^). 

On the other side, 



(3.13) 



(3.14) 



Q £ (XeU n ,XsUn) 



4(l-e 2 / 3 y) 1/2 \dy(xeu n )\ 2 + 



2 , W £ \XeU n \ 



dy 



(1 _ £ 2/3 y )l/2 

4/ (l-s 2 /* y ) 1/2 X?<dy + 8 [ {I - e 2 "y) l/2 X ! e Xe<u n dy 



+±£\l-s 2 / 3 y) 1 / 2 u' 2 dy + 4^ 



-2/3 



(1 - e 2/3 y ) 1/2 X £ V> 



+ 



-5/2 



-*/2 (1 " £ 2 / 3 y)V2 ^ T y £ -,/ 2< | yK ^ (1 - £2/3y)l/2 

The first two integrals in the right hand side of (|3.15p are 0(e 2//3 ), because u n E H 



dy. (3.15) 



\\x'Al~(j s )<s 2/3 



and 



max 



{(l- £ 2/3 y) l/2 :yeSuppX£ j^^3/2 



The fourth and last integrals in the right hand side of (|3.15|) are also 0(e 2 ^ 3 ), thanks to Lemma 
13. li From Lemma 13. 11 we also infer that 



/e r £ r+00 

(1 - e^yfl^dy = (1 + 0(e 2 ^)) / u' 2 dy = / u' 2 dy + 0(e 



_2/3-<5\ 



(3.16) 



^From Theorem Q] and from the decay properties of the function v n for n ^ 1 provided in (H n ), 
we deduce that v e = z^o + s 2 ^ 3 r £ , where r £ = O^oo (r)(1) and z^or e = 0z,°°(r)(1) as e — » 0. As a 
result, W e - W = 3(u 2 - v 2 ) G L°°(R), and 



\\W £ - W 1| L» OK) 



< p 2 /3 



Then, since Wo(y) = 0(y) as y — ► ±00, 



Wn 

(l- £ 2/3y)l/2 ° 

As a result, using once more Lemma 13,11 



< =-2/3-5 



L oo(_ £ -«/2 )e -i/2) 



(3.17) 



(3.18) 



-5/2 



W e ul 



-S/2 (1 - £ 2/3y)l/2 



f +00 



2/3-5 > 



(3.19) 
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Finally, we get from (|3.7p . f|3.14j) . (|3.15p . (|3.16p . (|3.19p and the estimates on the other term in 
the right hand side of A3. 15f> : 

R £ (v n ) = R(u n ) + 0(e 2 / 3 - 5 ) = ^ n + 0(e 2 / 3 - 6 ), (3.20) 

which completes the proof of (|3.5p and of its corollary (|3.6p . 



Step 2. Asymptotic behaviour of the eigenfunction u n . Property (i) n will be obtained 
as a consequence of (I3.6P and of the converse inequality 

/X n ^/2n + 0(£ 2/3 ~ <5 )- 

The proof of the latter inequality is delivered in Step 3 below. The proof uses the following 
properties of the eigenfunction u e n corresponding to the n th eigenvalue jl n of M £ . 

Lemma 3.2 There exists a constant C n > such that for every y £ J e and e > sufficiently 
small, 

K{v)\<C n e-M, (3.21) 

whereas 

( C n (\y\ + l)e~\y\ if y^O, 

C n (|y| + l)e-H+exp(-^) if < y < 



m)'(y)\ < 



(3.22) 



if £~ 2/3 < y < e -2/3 



Proof. In order to prove (|3.2ip . we come back to the eigenfunction 

(f E 2n (x) = u e n 



1-x 2 

of L £ + corresponding to the eigenvalue A|„ = fl n e 2 ^ 3 . Since 

ll^nll/fi = Q (^n)^n.) = Arall^ralle = P"n 

it follows from (|3.6p and Lemma 12.41 that for e sufficiently small, 



IbLIU-(M) = KIU»(J«) < C^/W n ^ C^Jn n + 0(e^~ 5 ) ^ c n , (3.23) 

where c n > is an e-independent constant. Since Wo(y) > \y\ as y — ► ±oo, we can fix a n large 
enough such that inf{Wo(y) : |y| ^ a n } > 4 + /i ra . Then, using (|3.17p and (13. 6p . we obtain, for 
x 2 < 1 — a n e 2 / 3 and for e small enough, 

(_ £ 2 5 2 + x 2 _ j + 3r? 2 _ ^2/3^) exp / I - ' J 



£ 



2/3 



= ^ 2/3 ("^ 2/3 " ^ + ^ " ft) «p (-^7?) 

;> e 2/3 (^_ 4 + inf{W / (y) : y ;> «„} - /i„ + 0{e 2 ' 3 ~ 5 )) exp ("^r) > 0. (3.24) 
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On the other side, <pf, n solves the differential equation 

(s 2 d 2 x + x 2 - 1 + 3r, 2 - e 2 / 3 faM n = 0. (3.25) 

Thus, 

{-e 2 d 2 x + x 2 - 1 + 3r? 2 - e 2/z fL £ n W n± > 0, |x| < (1 - a^ 2 / 3 ) 1 / 2 , (3.26) 

where 2 

^±(x) = c n exp ^a n - 1 g2/ 3 j ± </?l n (x). 

Moreover, from (|3.23|) . we get 

^ ± (±(l-a n£ 2 / 3 ) 1 /2)^ . 
As a result, since for e small enough, we also have like in (|3.24[) 



> e 2 / 3 (inf{W (j/) : y > a n } - » n + 0(e 2 ^- & )) > 0, (3.27) 
the maximum principle ensures that 

V£±(aO£0, |x| < (1 - a^ 2 / 3 ) 1 / 2 , 

which is equivalent to 

W £ 2n( x )\ < c n ex P f Q n ~ ^2/3 1 i M < (1 ~ O n £ 2/3 ) 1/2 . 

In terms of u £ n , it means that 

K(y)\ < c n e a "e-y , a n < y < e" 2 / 3 . (3.28) 
On the other side, for |x| ^ (1 + a^e 2 / 3 ) 1 / 2 and for e sufficiently small, we obtain like in (|3.24|) 

(- £ 2 8j + - 1 + 3,| - e 2 / 3 K) exp (-^J 1 ) 
= - ^ + Wc (l_Lf!) _ ft) exp (-^) 

^ 0. (3.29) 
Thus, exp f~ ^7r l i s a positive, continuous supersolution of 

(_ e 2 5 2 + x 2 _ 1 + 3r? 2 _ £ 2/3 -e ^ = q 
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in {x : \x\ > (1 + On^ 2 / 3 ) 1 / 2 }. From a slightly modified version of Corollary 2.8 in PQ, we deduce 
that 2 

\tp E 2n (x)\ ^ 2c n exp (l + a n - , \x\ > (1 + K + l)e 2/3 ) 1/2 . 

More precisely, the constant 2c n e^ an+l ^ above has been chosen in such a way that the inequality 
holds for \x\ = (1 + (a n + l)e 2//3 ) 1//2 , and the result in pj] ensures that then, the inequality holds 
for any x such that \x\ ^ (1 + (a n + l)^ 2 / 3 ) 1 / 2 . In terms of u^, it means that 



\u £ n (y)\ < 2c n e a "+V, y^-(a n + l). 



(3.30) 



Then, ([3311 follows from (13T23D . (I3T28|) and (l3T30|) . We next prove (1332D . From (133X1) and the 
differential equation M e -u^ = ji £ n u £ n , we infer that for every y £ J £ , 



d y {l-e^y) l,2 d y u £ n ) (y) 



C n e 



4(1 - {-2/3y)l/2 



^+(sup^ M+0(s 2 /^) 



(3.31) 



where we have also used (|3.6p and (|3.17p . The estimate (|3.22p in the case y < directly follows 
by integration of (|3.3ip between — oo and y: 



ffjm^ (l-e^y^^Jiy) 



dyil-e^yf^dyul) (s)ds 



< (|y| + l)e-K (3.32) 



— 2/3 —-2/3 

As for the case < y < £ —^ — , integration of (|3.3ip between y and — gives 



[ =--2/3 

(1 - e^yfl\u £ n )'{y) - 



(3.33) 



which provides thanks to the triangular inequality 



+ i) e -ifi + 



-2/3 



(3.34) 



Using basic integration, we also have 



e -2/3 ff -2/3 
«n(— g— ) - «^(— 4— ) 

.-2/3 / 



-2/3 » 



^ IK)'(S) _ e 2/3 s )l/2 



-2/3 (1 _ £ 2/3 s )l/2 



ds. (3.35) 



Since the last integral in the right hand side of (|3.35p is bounded from below by 
deduce from (|3^5|> . ([3331) and ([3311 that 



4 ' 



-2/3 
(^)'(^i-) 



exp(- 



-2/3 



-)■ 



(3.36) 
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—2/3 

Combining (j3.36j) and (|3.34p . we get (|3.22p in the case when < y < — • Finally, we consider 
the case when 6 2 /3 < y < e -2 / 3 . Integration of (|3.31|) between £ 2 /3 and y yields 



{u E n )'{y)\ < 



< 



< 



(1 - e 2 /3 y )i/2 

1 

(1 - e 2 / 3 y y/ 2 

£ -2/3 



-2/3 



+ 



(|s| + l)e-l s l 

2/3 (! _ £ 2/3 s )l/2 
e -2/3 



exp(- 



--2/3 c-2/3 

- 2/3 exp(-i— ) 



+ e 



-V3 (l_ e 2/3 s )l/2 



exp(- 



(3.37) 



(1 - £2/3y)l/2 ' 

where we have also used f|3.36j) . This completes the proof of (|3.22p and the proof of Lemma 13.21 



Step 3. Lower bound on jl £ n and proof of (i) n . In order to show that (i) n holds, we next 
prove the converse inequality 

Hn < fi £ n + 0(e 2 / 3 - 5 ), (3.38) 
which will be deduced from (|3,3p and 

n-l 

R(v £ n ) = ^ + 0(e 2/3 ~ 5 ), where v £ n = X eU £ n - ^ (XeU £ n , u k ) u k , (3.39) 

fc=i 

In order to prove (|3.39p . we proceed similarly as for the proof of (|3.5p . First, since (iii)k is assumed 
to be satisfied for k ^ n — 1, 

1-1 2 
Q(XeU £ n , Xe<) - E ^fc (XeU E n ,U k ) 

m) = ^ — 

llXe^nP - E (XeU £ n ,U k ) 2 
k=l 

Q{XeU £ n , XeUn) ~ Vn-1 (XeU £ n , U n _i) 2 + 0(g 2 / 3 ~ <5 ) 

||x £ ^H 2 -(^^n-i) 2 + 0(e 2 / 3 -' 5 ) 
Then, thanks to Lemma 13.21 and the normalization of 

g -2/3 

llX^nll 2 = / J 2/3 Xel^nl 2 ^ 
, ' -£ ~~2 — 

e |2 
nl 

a (l- £2/3^1/2-" 1 J £ - S<yK! ^l 



(3.40) 



(1 + 0(e 2 / 3 ^)) f ] % dy + f _ 2/3 xlKfdy. 



(^0(e 2 ^))J j Ai J% y)1/2 dy + 0(e 2 ^ 

l + 0{e 2 ^- s ). (3.41) 
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Similarly, using Lemma 13.21 and (|3.18p and proceeding as in (|3.15p , we get 



+00 



Q{XeU £ n ,XeU £ n ) = I (4\dy(x £ u £ n )\ + W \xsU £ n \ ) dy 

hoo r+00 



/+00 p+00 
X ' F 2 \u £ n \ 2 dy + 8 XeXe{u £ n )'u £ n dy 
-00 J —00 

+4(1 + <D(e 2 ^)) f \l-e^yfl\ul)fdy 

J-e-* 



+4 / 2/3 X 2 \(n £ n )'\ 2 dy 

f £ ~ S/2 f 

+ / W \u £ n \ 2 dy + W \xeU £ n \ 2 dy 
= 4(l + 0(e^- s )) [ {l-e^y) l/2 \{u £ n )'\ 2 dy 

JJe 

re- 5 / 2 w \~.e \1 

= Q £ (u £ n ,u £ n ) + 0(e 2 / 3 ~ s ) = fi £ n + 0(e 2 ^- s ). (3.42) 

In order to deduce (|3.39|) from (|3.40p . it remains to estimate the scalar product (xsU„, u n -\). 
Notice that in the case when n = 1, this term does not exists, and there is then nothing to do. 
From (iv) n _i, there exists c n _i S M such that 

\\ X eK-l ~ C^itin-iH < E 1 '^' 2 . (3.43) 

Then, by triangular inequality, and thanks to (|3.4ip for n replaced by n — 1, 

|| Cn _i| - IK \\Cn-lU n -l -Xe^-lH + |||X^n-lll ~ l| £ S 1/3 ~ 5/2 , (3.44) 

whereas 

\c n -l\ \ (XeU £ n ,U n -i)\ 

< \(XeU S n ,Cn-lU n -i -Xe«^_i)| + ( ( x\ ~ n_ J/3 y y/2 J "n» "n-1 ) + | <«n> «n-l> e | 



< e^' 2 , (3.45) 

where the first term in the right hand side of (|3.45p has been estimated thanks to the Cauchy- 
Schwarz inequality, (|3.43p and (|3.4ip . The second term has been estimated thanks to Lemma [3.21 
for u e n and for uf l _ 1 , and the last one is equal to 0. We deduce from (I3.44|) and (|3.45p that 

(Xeu £ n ,u n - l ) = 0(e 1 / 3 ~ s / 2 ). (3.46) 

Then, (|3^S9j) and (EOHl) follow from Ml . (I3T411 . (13^421) and TOT) . Property (i) n is a direct 
consequence of ()3.38|) and (|3.6p . 

Step 4. Proof of (ii) n and (iv) n . From the definition of v £ n in (|3.39p . it is clear that 

v £ n e Span(ni, • • • ^Un-x)^ . 
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Thus, v e n can be decomposed as 

v £ n = c e n u n + w £ n) where c^£l and w £ n S Span(m, • • • , it,,) 1 . (3.47) 
From (|3.39p and (i) n , we have 

^ + o( e / ) = Mn + o( £ / (c , )2 + | K ||2 > + w ■ 

It follows that 

(/i„, +1 - Mn )ll<l| 2 <e 2/3 "1^H 2 - (3-48) 
Thanks to the definition of v e n in (|3.39p . property (iii)k for k ^ n — 1 as well as (|3,46p . 



~e || < .1/3-5/2 



(3.49) 



On the other side, (|3.4ip ensures that Hxe^nll -> 1 as e -> 0. As a result, \v^\ — > 1 as e — ► 0, 
and (|3.48p implies 

IKII <e 1/3 - 5/2 . (3.50) 
Moreover, from Lemmas 13.11 and 13.21 we infer that for any k ^ 1, 

(XeU e n ,u k ) = [ XeU £ n u k dy+ f _ 2/3 xSn u h dy 



= (1 + O(e^S)) {l ^f; )1/2 dy + O(e^). (3.51) 

^From (|3.50p and (|3.5ip we deduce in particular that for every k n + 1, 

0(e 2 / 3 ) + (x^k, <) £ (1 + 0(e 2 / 3 - 5 )) = ( X eu £ n , u k ) = (v* n ,u k ) = «, u fc ) = ©(e 1 / 3 " 5 / 2 ), (3.52) 
which proves (ii) n - Then, (iv) n is a consequence of the triangular inequality, (|3.49p and (|3.50j) : 

\\ Xe U £ n ~ C £ n U n \\ ^ \\ X eU £ n ~ V £ J + ||<|| < e 1 ^ 2 . 

Step 5. Proof of (iii) n and (v) n . Like in (|3.47p . we decompose v s n _ x as 

v n-l = g n-l<-l + Wn-li where d n-i G R and G Span(ttf , • • • ,u*_l) ±e - 

^From (|3.5p for n replaced by n — 1 and (i) n -i, we have 



(^-i) 2 + IK-ill 2 

(^n-l) 2 An-l + H^n-llleAn 



(*£-i) 2 + K_ill§ 
Using (i) n and (i) n -i, it follows that 

( Mn - + o^ 2 / 3 - 5 ))!!^!!! 2 = (ft - ^-OIK-xll 2 < e^Wv^wl (3.53) 
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Thanks to the definition of vf l _ 1 given by (|3.5p and property (ii)k for k ^ n — 2, 



\ v n-l ~ Xe u n- I !- 



i|| £ < £ 1/3 - 5/2 . (3.54) 



Thanks to (|3,14p for n replaced by n— 1, HXe^n-llle — > 1 as e — > 0, thus 1 1 f ^ 1 1 1 as e -> 0. 

As a result, we deduce from (]3.53j) that 

IK-il| £ <^ 1/3 - 5/2 . (3-55) 

Then, for every k ^ n, we get 

(Xeul,un-i) (1 + 0{e 2 ^- 5 )) = (xeu n -i,u £ k )e = «-!,&%) e = «-!,**) 6 = 0(e 1/3 ^ /2 ),(3.56) 
using similar arguments as in the derivation of f)3.52p . Moreover, 



(Xeu n -i,u e k ) £ = (1 + 0(e 2/3 5 )) I {x e u%,Un-i) - / 2/3 XeU n ~iu%dy 



r XsUn-lUj 
^-^l^/i (l- £ 2/3 y) l/2 » 

= (l + 0( e 2 / 3 - 5 ))( Xe n|,n n _ 1 ) + 0(e 2 / 3 ), (3.57) 

where the two integrals in the right hand side of (I3.57|) have been estimated thanks to the Cauchy- 
Schwarz inequality, Lemma |3. II and the normalization condition ||u||| e = 1. The combination of 
(I3.56P and (13.571) completes the proof of (iii) n . Then, (v) n follows from the triangular inequality, 
(|53aj and (I335|) : 

||X £ n n _i - ^_i^_i|| s ^ \\XeUn-l ~ «ft_i|| e + IK-llle < £ 1/3 ~ 5/2 . 

It completes the proof of (G n ), and therefore the proof of Theorem [2 ■ 

4 Semi-classical limit for eigenvalues of L £ + 

We list here formal results of the semi-classical theory that describe the distribution of eigenvalues 
of L e + . We will show that the standard Bohr-Sommerfeld quantization rule does not give the 
correct asymptotic behavior of the eigenvalues of L £ + as e — ► because the potential V e (x) depends 
on e. Nevertheless, the Bohr-Sommerfeld quantization rule gives the correct scaling 0(e 2//3 ) in 
agreement with the asymptotic limit (|3.2p in Theorem [2j 

Eigenvalue problem for operator L e , can be rewritten in the form 

(-d 2 x + e~ 2 V e {x)) u{x) =s- 2 Xu(x), (4.1) 

By properties of rj e following from Theorem [H the potential V e (x) has the properties 

• V £ (x) G C°°(R) for any small e > 0, 

• lim!4(a;) = Vq(x), where Vq £ C(M) is given by 



V Q (x) 



2(1 -x 2 ), \x\^l, 
x 2 - 1, \x\ ^ 1, 
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• V e (x) takes its absolute minimum at ±a e for any small e ^ and a £ — ► 1 as e — > 0, 

• V^(x) — ► +oo as |x| — ► oo for any small e 0. 

If F e (x) is replaced by Vo(x), the eigenvalue problem (|4.ip takes a simplified form 

+ e~ 2 V (x)) u(x) = s~ 2 \u(x), xeR, (4.2) 

which describes the eigenvalues of the operator L £ + mentioned in section 1. As it is well-known 
(see a recent review in [BDS| ). the eigenvalues of the Schrodinger operator —d 2 + e~ 2 V(x), with 
a smooth, e- independent double well potential V(x), are twice degenerate in the semi-classical 
limit e — > 0. Namely, the eigenvalues are grouped by pairs. In each pair, the two eigenvalues 
are exponentially close one from another as e — * 0. The asymptotic distribution of these pairs of 
eigenvalues is determined by the Bohr-Sommerfeld quantization rule. 

Let us try to apply the Bohr-Sommerfeld quantization rule to the eigenvalue problems (|4.ip 
and (|4.2p for the operators L e + and L e + , in spite of the fact that this rule was proved rigorously by 
Fedoryuk [Fed] only for a class of e-independent, analytic potentials. Since neither (|4.ip nor (|4.2|) 
satisfies assumptions of the main theorem in [Fed] . this application is purely formal. According to 
the standard Bohr-Sommerfeld rule, the consequent eigenvalues A| n _ 1 and Af n of the Schrodinger 
equation (14. ip with the double- well potential V e {x) would be given asymptotically by 



r*%W / i\ 

/ i/A - V £ (x)dx ~ £7r I n , as e — > 0, for fixed n ^ 1, (4.3) 

JxtW V 2/ 



where x^(A) are the roots of V e (x) = A on R + , such that < x £ __(X) < 1 < x+(A) < oo. Let us 
use the scaling 

y = ^AT> V £ {x)=e 2 ^W £ {y), \ = e 2 '^, (4.4) 

where W £ (y) = 3v 2 (y) — y and /i is a new eigenvalue. The Bohr-Sommerfeld rule is rewritten in 
an equivalent form by 

rv+(p) y/n - WAy) 

I dy ~ 7r(2n — 1), as e — > 0, for fixed n ^ 1, 

V(m) Vl-e 2/3 y 

where y±{n) are the roots of W £ (y) = /i on I, such that — oo < yl(/u) < < y £ + {^) < oo. Taking 
the limit e — > for a fixed n ^ 1, we obtain 



/•j/+(a0 

/ v / / i -^o(y)dy~vr(2n-l), for fixed n > 1, (4.5) 



where Wo(y) = 3fg(y) — y and y±(/x) are the roots of Wo(y) = M on K. The new expression is the 
Bohr-Sommerfeld quantization rule for the Schrodinger operator Mq = — 4<9 2 + Wo and it is only 
valid for large n> 1. Therefore, the Bohr-Sommerfeld quantization rule (14. 3p does not recover 
the statement of Theorem [2] correctly. Meantime, it still implies that the eigenvalues A| n _i and 
A| n for a fixed n ^ 1 are scaled as C(e 2 / 3 ) as e — > 0. The discrepancy of the Bohr-Sommerfeld 
rule is explained by the fact that the smooth potential V £ (x) in the eigenvalue problem (|4.ip 
depends on e. 
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Note that the limit e — > can be computed exactly for the simplified eigenvalue problem (|4.2p 
thanks to the scaling transformation (14. 4h , In this case, the limiting formula (|4.5h holds with 
Wo(y) replaced by 2y for y ^ and — y for y ^ 0, so that y-(fJ-) = —fj, and y+(n) = fJ>/2. In other 
words, 

/•o /-m/2 

/ + + / vA* ~~ ~ 7r(2n — 1), for fixed n ^ 1, 

and the computations of integrals gives /in ~ (7r(2n — 1)) 2//3 , in agreement with the behavior 
0(n 2 / 3 ) of eigenvalues of the Schrodinger operator with a linearly growing potential as |y| — > oo 
|Suj . Therefore, the Bohr-Sommerfeld quantization rule suggests that the eigenvalues {A^} n ^i of 
the simplified operator L £ + considered in our previous work [GPJ satisfy the asymptotic limit 

lim =lim% = (vr(2n - 1)) 2/3 , for fixed n > 1. (4.6) 

e|0 £ 2 / 3 e|0 £ 2 / 3 

However, the justification of the asymptotic limit (14.60 cannot rely on the work of Fedoryuk [Fed] 
because the e-independent potential Vq(x) in the simplified eigenvalue problem (|4.2p is continuous 
but not C 1 on R. 



5 Proof of Lemma 12.11 

Let a > 1 be like in the assumption of the lemma, and A = \\x a /||l°°(k + ) < °°- We first prove 
(12. lip by contradiction. We proceed as follows. We suppose that (12. lip is not true. Namely, we 
make the assumption 

(G a ) <p(x) + 0(x-( Q+1 )), x a f £ L°°(A + , +oo) 

If a > 2, we prove that (G a ) implies (G a _2), such that after a finite number of steps, (G Q ) 
implies (Ga) for some a G (0,2]. On the other side, we show that for < a ^ 2, (G a ) yields to 
a contradiction. 

If (12. lip is not true, then, up to a change of / and ip into — / and — cp, there exists a sequence 
(x n )n^n (where uq > A), such that x n j oo, x n ^ A + and 

x%W(x n )ip(x n ) > n. 

Then, 

x"ip"(x n ) = x"W(x n )ip(x n ) - x"f(x) ^ x%W(x n )ip(x n ) - A > n- A. 
For n ^ no > A, we define 

2/n = sup{y > x n ,Vx G (x n , y), x a W(x)ip(x) - A > (n - A)/2}. 
By continuity of W and <p, for every n ^ no, either y n = +oo or 

. . Th ~\~ A. 

We distinguish the two following cases: 

A) There exists n\ ^ no such that y ni = +oo 

B) For every n ^ no, y n < +oo. 

In case B), extracting a subsequence of (x n ) n j> no if necessary, one can assume that 
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For n ^ no + 1, we define 

x„ = inf{y < x n ,Vx G (y, x n ), x a W{x)<p{x) - A > 3(n- A)/4}. 

Since y n -i < x n and 

yS-iWWlMlto-i) - A = (n - 1 - A)/2 < 3(n - A)/4, 

we deduce x n > y n -i > — oo. Moreover, by continuity, ip(x n ) = (3n + A)/(4x%W(x n )), and 
(f(x) > (3n + A)/(4x a W(x)) for x > x n , x close to x n . Therefore 

,...371 + 4(1 / 1 \ ^ 3n + A/ \\W'\\ L °°\ 1 



4 ds V^N/l^ 4 V C + J x« +l W(x n ) 

> —C^n n , 

for some Ci > 0. By definition of y n and x n , for every x G (x n ,y n ), 

x <p (x) ^ ■ 

Thus, 

(x) > p (x n ) + — — / — 

1 77 _ A r x 1 
> -Cm—-, + / — dy =: GJx). (5.2) 

Notice that G n (x n ) < 0, whereas 

+oo if a ^ 1 



Gn(+00) 

where for a > 1, 



# n if a > 1, 



~ — , x~ a -i > as n ^ +oo. 



n - A 
2(a-l)x£ 

As a result, for n sufficiently large, since G n is increasing on (x n , +oo), G n vanishes exactly once 
on that interval. Moreover, this unique zero z n of G n is defined by 



1 , 2dn 1 
— dy 



y a y n-Ax% +1 W(x n ) 
thus 

= X n + O ^^2^ • 

By integration of (|5.2p . we infer that for x G (x n ,y n ), 

ip{x) ^ y?(x n ) + / G n (y)dy 

> </?(x n ) + / G n (y)dy 



3n + A 1 C 2 n 



4 5%W(x n ) x£ +3 W(x n )' 
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for some constant C2 > 0. Therefore, for n large enough, for every x G (x n ,y n ), since W is 
increasing on (A+,+00), 

5n 1 5n 1 

For n sufficiently large, 5n/8 > (n + A)/2, and it provides a contradiction with (|5.ip . which 
means that case B) can not happen. In case A), for every x x ni , 

iVWM»i-^)/2>0. (5.3) 

Therefore <p'(x) | as x | 00, otherwise 92 would not be in L 2 (M). Thus, for every x ^ x ni , 
(p'(x) ^ 0, and therefore <p{x) [ as x \ 00. If < a ^ 1, (|5.3p provides a contradiction with the 
fact that <p'(x) as x -> 00. If a > 1, integration of (|5.3p between a; and +00 yields 

^P^z 1 -. (5 - 4) 
z(a — 1) 



This is a contradiction with (p(x) — > 0, if 1 < a ^ 2. Finally, if a > 2, by integration of (|5,4p . 

fix) ^ / y dy = -x . 

2 {o>-l)J x 2(a-l)(a-2) 

Thus, 

<p(x) + olx-^- 2 ^ 1 



Since the assumption x a f G L°°(A + , +00) implies x a ~ 2 f G L°°(A + , +00), we have proved that 
(G Q ) implies (G5) if a > 2. The proof of (|2.1ip is completed by induction. Then, since (p" = 
Wip — /, we deduce 

<p"(x) = 0{x~ a ). (5.5) 

We next prove that 

^ (X) = | (1) if0<a<l. (5 - 6) 

By integration of (|5.5p . if a > 1, 92' (x) has a limit as x — » +00. This limit can only be 0, because 
92 G L 2 . (|5.6p is then obtained by integration of (|5.5j) between x an +00. If a ^ 1, (|5.6p is a 
consequence of the fact that tp(x) — > and tp"(x) ->0asi^ +00. 
Let x G C°°(R) be such that 

/ \ _ / if x < 1 
~ \ 1 if x ^ 2 

For m G N, let 9? m , / m G C°°(R) be the functions defined by 

<p m {x) = X (x)x-^ m+1 ) 

and 

frn(x) = -^(x) + W{x)if m {x). 
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^From now on, we assume that / and W have asymptotic series (|2.12p as x — > +00, so that 

+00 

f m (x) « x- (a+7m) yv fc x- 7fe + (a + 7m + l)(a + 7m + 2)x- (a+7(m+3/7)) . 

.7? — ► 4-rsn ' ^ 



k=0 

+00 

k=0 



where = \i k ^ 3/7 and u 3 / 7 = i> 3/ / 7 + (a + 7m + l)(a + 7m + 2). Notice also that the 
assumption W(x) ^ C + x implies vq ^ C + > 0. As a result, there exists coefficients {c ni ) m &% 
such that for every M ^ 0, 

M 

m=0 

where gM(%) = 0(x _a_7 ^ /+1 )) as x — > +00. Then, 

Af 
m=0 

where V>M = (— 9% + W) _1 yA/- Thanks to (|2.1ip . (|5.6p and (|5.5p . for M large enough, ^m{x) = 
Q^-a-^M+l)-!^ ^> M (x) = 0( x -«-7(M+i)+i) and ^(x) = 0{x~ a -^ M+ ^). Since this is true 

for arbitrarily large values of M, then (|2.13p and (|2.14p follow. 



6 Proof of Lemma 12.21 

By Proposition 12.11 we know that v$ is a strictly increasing function on R, with asymptotics at 
±00 given by (|2.9p and (|2.10p . Moreover, vq has a unique inflection point. From the behaviour 
of vo(y) as y — » ±00, we infer that Wo(y) = 3z^o(y) 2 — y — > +00 as y — ► ±00. We are going to 
prove that the global minimum of Wq is actually strictly positive. We argue by contradiction. If 
it is not the case, we can define 

Vi = inf{y > 0, u (y) = \/y/3~}, 

where we recall that Wo(y) > if y ^ 0. By continuity, v 0(2/1) = \J yi/3. We also denote the 
unique inflection point of uq by yo- Since f > solves (|2.5p . yo > is the unique solution of the 
equation ^o(yo) = \/y0) and ^ '(y) > if 2/ < 2/o> whereas v'^iy) < if y > yo. Notice that since 
fo(0) > and fo(yi) = \/yi/3 < Jyi, we have necessarily < yo < 2/1- Moreover, since uq is 
strictly increasing, we have y^yo = ^o(yo) < ^o(yi) = y ?/i/3, and therefore < 3yo < yi- 

First step: upper bound on yi. For y > 0, we introduce the function z(y) = VQ(y)/y/y 
and rewrite (|2.5p in terms of z(y) as 

Ay) + -'W = ^U) 2 -i + 4)- 
y 4 v y J 

Since z(y) — > +00 as y — > + and z(y) — ► 1 as y — > +00 with z(y) < 1 for y large enough (because 
for y > yo, v'o(y) < and therefore z-'o(y) < y/y), we deduce that z(y) admits a global minimum 
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at y = y m > 0, where 



< z"{y m ) 



IJrnZijJr, 



z(y m ) 2 - 1 + 



1 



4 \ v° 

\ an 

The assumption of non-positivity of Wq implies that z(y m ) ^ l/v3- Thus, 



i/m ° 



As a result, since u (y m ) ^ y/y m /3, 



3yo < 2/1 ^ 2/m ^ 



1/3 



(6.1) 



Second step: upper bound on ^(yo). Since is increasing on R and fo(y) 2 — y > if 
y < j/O) w e deduce, for every y < yo, 



I/O 



^o(t) 



[v Q {tf-t)dt 



(M m f-t) dt = ^(y -yf. 



(6.2) 



By integration, it follows that for y < yo, 



I/O 



v Q (t)dt 



< Vyo- ^o(yo)(yo - 2/) + -^rG/o - y) 3 - 



(6.3) 



The right hand side reaches its minimum (for y < yo) at y = y p , where y p < yo is defined by 
{yo ~ y P ) 2 = 8vo(yo)/y/yb~, and (HDI &ty = y p yields 



My?) < \/yo ■ 



4V2^ (y ) 3/2 



1/4 
2/o 



Since t'o > 0, the right hand side has to be strictly positive. Therefore 

1/3 



f9\' 



(6.4) 



Third step: upper bound on v' (yi). On the one side, notice that for y > yo, v'o(y) < 0, 
and therefore v (yi) ^ u' (yo)- On the other side, if y < y\, fo(y) 2 > y/3, and fo(yi) 2 = 2/i/3, 
thus 



^0(2/1) < 



d /y 



<iy V 3 



v=v\ 



1 

VSyT 
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As a result, thanks to (16.41) and (16.1 



**> * -((Iff 

Fourth step: upper bound on v' Q (y) for y > y^. For 6 G (0,2/3) to be fixed later, we 
define 

y 2 (6) = sup{y>j/i,Vt€ (yuy),^? < (1 — 5)t} 
(notice that v Q (y{) 2 = y 1 /3 < (1 - <%i). Then, for every y G (yi,y 2 (6)), 

= + f y ^(Mt) 2 -t)dt 



< v' (yi)+ [ yu -^(-st)dt 

Jyi 4 

< ^(yi) - ^pW - y?). (6-6) 

Fifth step: bound from below on y?(6)- For 6 G (0,2/3), we introduce the function h$ 
defined for y > y\ by 

^From (|6.5p and since 5 < 2/3, we infer h$(yi) < 0. Thus, if we define 

y 3 (<5) := sup{y > yi,Vt G (y 1 ,y),h 8 (t) < 0}, 
we deduce from (I6.6P that for y G (yi, min(?/2(<5), 2/3(^))) ; 

*>o(y) - Vi-<tyy = ^o(yi) - Vl - 6^/yl + [ V h s (t)dt < 0, (6.7) 

which implies that 

2/3(5) < ^(ff). (6.8) 

Sixth step: y 3 = +00. We shall see next that for an appropriate choice of 6, 1/3(6) = 
+00, which implies that y 2 (6) = +00 thanks to (|6.8p . This provides a contradiction with the 
assumption of non positivity of Wo, since vo(y) ~ yjy as y — > +00. An elementary calculation 
shows that /i^ reaches its maximum (for y > yi) at 

( VWT=6 \ 2/5 
y = yM '- = \-^yl6-) >yi ' 

where the inequality comes from (|6.ip and from the fact that 5 < 2/3. From (16. 5p . we obtain 

M»«) < ™ (7 » V"^ ' ) + *V - ^^g. (6.9) 
V ; \\32J v/3 2\/%I/ 8^3 8 -3Vio v ; 

For 5 = 1/3, elementary calculations show that the right hand side in (|6.9p is strictly negative for 
any y\ G (0, (3/2) 1 / 3 ), which implies that 2/3(1/3) = +00 and completes the proof of the lemma. 
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7 Proof of Lemma 12.51 

We denote 

U e {z) = W {e- 2/3 - e 2/3 \z\ 2 ), z e R d . 

We are going to show that there exists a constant C > such that for e > sufficiently small, 
for every ball B cR d , 

max U e (z) < — , / U e {z)dz. (7.1) 
zeB \B\ J B 

According to Theorem 0.3 in [Sh], Lemma 12.51 follows. First, we notice that, thanks to Lemma 
[2]2]and ([2. 150 . there exist C\, C 2 > such that for every y E R, 

Ci(l + |y|) < W (y) < C 2 (l + |y|). (7.2) 

Given zq E R d and r > 0, as z describes 5(zo,r), |z| describes the interval [|zo| — r, |zo| + r] if 
|zo| ^ r and the interval [0, \zq\ + r] if |zq| ^ r. Since the function 

/(s) = | e -2/3_ £ 2/3 s 2| ! seR+ 

is decreasing on [0,e -2 / 3 ] and increasing on [e _2//3 , +oo), we infer that max{/(|z|),z E B(zQ,r)} 
can only take the three different values depending on zq and r: either 

max /(|z|) =e 2/3 (|z | +r) 2 -e~ 2/3 and \z \ + r ^ e _2/3 (easel), 

or 

max /(|z|) =e" 2/3 -e 2/3 (|z | -r) 2 and < |z | - r < e~ 2/3 (case 2), 

z£B(2 ,r) 

or 

max /(|z|) = e~ 2 / 3 and |zo| — r ^ (case 3). 

ze-B(z ,r-) 

We are next going to prove (|7.1|) in each of these 3 cases. 

Case 1. We first show that for every zo, r like in case 1, we have 

N + -j= > e" 2/3 . (7.3) 

Under the extra assumption 

£ 2/3 ( | Zo | +r) 2_ e -2/3^ ff -2/3 ; (7 _ 4) 

(|7.3p clearly holds. On the other side, if (|7.4p is not true, then |zo| — r > since otherwise, 
E [|zo| — r, |zo| +r] and max{/(|z|),z E -B(zo,r)} /(0) = e -2 / 3 , contradicting the assumption 
that we are in case 1. Then, we also have 

£ 2 / 3 (|z | + r) 2 - e~ 2 l 3 = /(|z | + r) > /(|z | - r) > e" 2 / 3 - e 2 / 3 (|z | - r) 2 , 

which can be rewritten as 

\z \ 2 + r 2 > e^' 3 . 
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Since r < \zo\, we deduce 

£" 4/3 < \z \ 2 + r 2 \z \ 2 + r - + V2r\z \ = (\z \ + -^=\ , 
which means that (H3|) also holds if (|T1| is not true. Let a = y/3/2 ^ 1/2 + y/2/4. Then, 
(, 2 / 3 (N + arf - e~^) - \ (e 2 l\\z,\ + r) 2 - e~ 2 / 3 ) - \ (>( N + -^=) 2 - 
= 2e 2 / 3 \z \r(a-±-^=) + e 2 ^r 2 (a 2 -1/2-1/4). (7.5) 
We deduce from (|7.3p and (|7.5|) that for every z S B(z ,r) such that |z| > |z | + ar, 

/(|*|) = e 2/3 |z|2 _ e -2/3^ e 2/3 (|zo|+ar) 2_ e -2/3 

> J (e 2/3 (|*b| + r? - e- 2 ^) = \ max f(\z\). (7.6) 

2 V / 2zeB(z,),r) 



Then, we conclude thanks to (17.21) and (17.61) that 



■ v. / > T^Td / (1 + /(N))^ 

> rJo / (i + /(N))^ 



> \Sh /r fl + /(N))^ 

J '{zeB(z ,r): z-^\z \+ar)} 

^ max (l + /(|z|)) 

2\M d \ zeB(z , r y Jyi UJ 

where v a denotes the volume of {z G £>(0, 1) : Z\ ^ q}. 
Case 2. The assumption that we are in case 2 implies 

£ -2/3 _ £ 2/3 ( | Zo | _ r) 2 = /(|zq| _ r) ^ /(|zq| + r) ^ £ 2/3 ( | Zo | + r) 2 _ £ -2/ 3j 

and thus 

ko| 2 ^No| 2 + r 2 < £ - 4 / 3 . 

It follows that 

^-2/3 _ £ 2/3 (ko| _ r/2) 2^ _ 1 ( £ -2/3 _ £ 2/3 (N _ r)2 ) = 1 ^-2/3 _ ffi^ + 

We deduce that for every z £ B(zo,r) such that \z\ ^ |zo| ~~ r /2, 

/(W)>^/(N-r). (7.8) 



4/3 2 
— > 0. 
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Then, we show that this last estimates holds as soon as z £ B(zQ,r) and z ■ zqI\zq\ ^ |zo| — 7r/8. 
Indeed, under this assumption, Pythagoras' theorem ensures that 

I 1 2 / z V . I 12^/ x Z N ~ 

/ . 7r \ 2 9 / 7r x 2 
< ko --5- +r 2 - ' 



r\ 2 3r / r 

^1-2) -T (N - r) n N -2 

Then, we conclude similarly as in case 1, thanks to (|7.2p and (|7.8p 
1 1 U £ (z)dz > / {l + f(\z\))dz 



B{z a ,r)\ J B (z ,r) " \^ d \r d Jb(z ,t) 



> J&sL Al + H\*\))dz 



[zeB(z ,r): z-^^\z \~7r/sj 



ClVj/8 

2\B d \ ziB(z , r y 



Case 3. First, we notice that the assumption that we are in case 3 yields 

e- 2/3 > /(M + r) > e 2 '\\z \ + rf - e^, 

which gives 



Thus, since \zq\ ^ r, we get 



If the extra assumption 



\zq\ +r < V2e~ 2/3 . (7.10) 



\zo\ < (7.11) 



r > 5\z \/4 (7.12) 



holds, then ()7.1ip and the triangular inequality give B(0,r/5) C B(zo,r). Moreover, if z E 
5(0, r/5), then we get from (|7.10p 

e -2/3 _ £ 2/3| z |2 ^ 23e" 2/3 /25. (7.13) 
Then, we conclude similarly as in cases 1 and 2: 

, , 1 / U £ (z)dz > [ (l + 23e~ 2/3 /25 N ) dz 

\B(z ,r)\J B{zo>r) £K J ' \M d \r d J mr/5) \ 1 J 

> ^r^r ^M*)- ( 7 - 14 ) 

ZO • 3 a U2 zt=B(z ,r) 
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As for the last case when (|7.12p is not true, we have then 

z G B(z Q ,r) : z ■ ^ ^ -J \ C B(z , r) n 5(0, |z |). 
M 40 J 

Indeed, using also \zq\ ^ r ^ we have then 

I |2 Z' Z V , I |2 A I z V I |2 I |2 i ol I Z 

^ 2 I 1 2 i ? r z - I |2 / 71 r\ 

r ~ ' °' ^o - ' ' ^ 

On the other side, for z G -B(0, |zo|)) thanks to (|7.1ip . we have 

/(|^|) =£ - 2 /3_ e 2/3| z |2 >£ -2/3 /2 _ 

Then, we conclude similarly as in the previous cases: 



[ U £ (z)dz > -^-j [ 



lp , „ , - - - .p.,,, , » (l + /(|z|))cfe 



^ iWR ft (i + f{\z\))dz 



^From ((EH), (fT5]) . (|7T4l) and ([77T6l) . we infer that (JUJ) holds, with 

r= . / ClUq gl^7/8 23Ci Ci ^33/40 

~ mm V2|B d |C2'2|]B d |C2'25-3 d C2' 2|B d |C 2 ^ ' 
which completes the Proof of the lemma. 
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